function [Wp, Wm, Qp, Qm] = WQ_n(~, x, z0, zc, n, p)
%WQ_N Summary of this function goes here
%   Detailed explanation goes here
    jn = besselj(n, x);
    jnp= besselj(n+p, x);
    jn_p = besselj(n-p, x);
    zn = z0+n*zc;
    lrtz = 1.0/(1.0+ zn*zn);
    lrtz_diff = zn*lrtz;
    
    Wp = jn.*jnp*lrtz_diff;
    Wm = jn.*jn_p*lrtz_diff;
    Qp = jn.*jnp*lrtz;
    Qm = jn.*jn_p*lrtz;
end

